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ABSTRACT 

The discovery of the eccentric binary and millisecond pulsar PSR 
J1903+03273 has raised interesting questions about the formation mechanisms of 
this peculiar system. Here we present a born-fast scenario for PSR J1903+03273. 
We assume that during the supernova (SN) explosion that produced the pulsar, 
a fallback disk was formed around and accreted onto the newborn neutron star. 
Mass accretion could accelerate the neutron star's spin to milliseconds, and de- 
crease its magnetic field to ~ 10 8 — 10 9 G, provided that there was sufficient mass 
(~ O.1M ) in the fallback disk. The neutron star became a millisecond pulsar af- 
ter mass accretion terminated. In the meanwhile the binary orbit has kept to be 
eccentric (due to the SN explosion) for ~ 10 9 yr. We have performed population 
synthesis calculations of the evolutions of neutron stars with a fallback disk, and 
found that there might be tens to hundreds of PSR J1903+03273-like systems in 
the Galaxy. This scenario also suggests that some fraction of isolated millisecond 
pulsars in the Galactic disk could be formed through the same channel. 

Subject headings: accretion, accretion disks — pulsars: general — pulsars: indi- 
vidual (PSR J1903+0327) — stars: neutron 

1. Introduction 

There are now around 100 binary and millisecond pulsars (BMPSRs) known in the 
Galaxy. Most of them are characterized by short spin periods (P < 15 ms), weak magnetic 
fields (B ~ 10 8 — 10 9 G), and circular orbits with companions of mass ~ 0.15M© — O.45M . 
In the recycling scenario, these systems are thought to be the descendants of low-mass X- 
ray binaries (Bhattacharya & van den Heuvel 1991; Tauris & van den Heuvel 2006). The 



institute of Astrophysics, Huazhong Normal University, Wuhan 430079, China 
2 Department of Astronomy, Nanjing University Nanjing 210093, China 



-2- 



evolutionary path is briefly described as follows. A high-field (B ~ 10 12 — 10 13 G), rapidly 
rotating neutron star (NS) is born in a binary with a low-mass (~ 1M ) main- sequence 
(MS) companion star. During the supernova (SN) that produces the NS, mass loss and a 
kick imparted on the NS cause the orbit to be eccentric. The NS spins down as a radio pulsar 
for ~ 10 6 — 10 7 yr until passing by the so-called "death line" in the magnetic field - spin 
period (B — P) diagram (IRuderman fc Sutherland! Il975l ). When the companion evolves to 
overflow its Roche lobe, mass transfer occurs by way of an accretion disk, and tidal friction 
serves to circularize the orbit. Mass accretion onto the NS gives rise to X-ray emission, 
induces magnetic field decay (the mechanisms for the field decay induced by accretion are, 
however, not well understood), and accelerates the NS up to short (millisecond) period. 
When the companion loses almost all of its envelope and mass transfer ceases, the endpoint 
of the evolution is a circular binary containing an NS visible as a low-field, millisecond radio 
pulsar, and a He or CO white dwarf (WD), the remaining core of the companion. A fast- 
growing populat ion of eccentric (e > 0.1) BMP SRs have been recently revealed in globular 
clusters (GCs) ( iFreire. Ransom fc Guptal 120071 ). These high eccentricities are likely to be 
attributed to the dynamical interactions in the central regions of GCs. 



Most recently I Champion et al.l (120081 ) reported the discovery of an eccentric BMPSR 
PSR J1903+0327 in the Galactic plane. With a spin period of 2.15 ms, the pulsar lies 
in an eccentric (e = 0.44), 95-day orbit around a ~ 1M© companion. The mass of the 
pulsar was estimated to be 1.74 ± 0.04 M , about 30% larger than that of other binary 
NSs in the Galactic disk. Infrared observati ons identified a possible MS companion star . 



Preferred formation scenarios may include (cf. Ivan den Heuvel 



20081 : IChampion et al]l2008h : 



(1) the pulsar was recycled in a GC, then the original donor star was replaced by the present 
MS companion via one or more exchange interactions, and the binary was displaced from 
the GC due to either ejection or the disruption of the GC; (2) the pulsar is part of a 
primordial hierarchical triple system, recycled by accretion from the progenitor of a massive 
(~ 0.9 — 1.1 Mq) WD which is seen in the timing measurement, while the detected infrared 
counterpart, the third star, is in a much wider and highly inclined orbit around the inner 
binary. The eccentricity of the inner binary is caused by the perturbation of the outer 
MS star; (3) in an alternative triple-star model, the WD in the inner binary system was 
evaporated by the pulsar's energy flux, or coalesced with it, such that the present binary 
remained. 

Since there is no strong evidence for the GC or triple star origin of PSR J1903+0327, in 
this paper, we seek an alternati v e, bor n-fast scenario for PSR J1903+0327 which was already 
mentioned by IChampion et al.l (120081 ) . These authors argued that the pulsar was not likely 
to be formed spinning rapidly at the time of SN with a small magnetic field, due to the 
following reasons: (1) there are no pulsars like PSR J1903+0327 in any of the > 50 young 
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supernova remnants in which an NS has been inferred or detected directly; (2) the 18 isolated 
MPSRs detected in the Galactic disk have spin distributions, space velocities, and energetics 
indistinguishable from those of recycled BMPRSs and their space velocities and scale heights 
do not match those of non-recycled pulsars; (3) magnetic fields in young pulsars is unlikely to 
be less than 10 10 G. However, as we shown below, if experienced accretion from a SN fallback 
disk, the newborn pulsar could be accelerated to a period of several milliseconds, together 
with the magnetic field decayed to ~ 10 8 G from an initial value, say ~ 10* 2 G. When the 
disk becomes neutral after ~ 10 3 — 10 4 yr (IMenou. Perna fe Hernquistll200ll ) and accretion 
ceases, the NS becomes an MPSR. The orbit, if not being disr upted by the SN explosion , 
will remain eccentric for > 10 10 yr before it is circularized (cf. iHurley. Tout fe Polsl 120021 ) . 
The binary consists of a recycled pulsar and a MS companion in an eccentric orbit, just like 
PSR J1903+0327. Hence this scenario could produce a population of eccentric BMPSRs as 
well as isolated MPSRs. 



2. Spin Evolution of An NS with A Fallback Disk 



apse of its pr o genitor, a small 



Chevalier 1989 



Following the formation of an NS through the core col 
amount of mass could fall back on to the compact object fjColgatd 11971 
Lin. Woosley fe Bodenheimerlll99ll ). Some of the fallback material may carry sufficient an- 
gular momentum and to form an accretion disk of mass M^ around the NS. The initial 
transient phase la sts for a local viscous timescale t n ~ 6.6 x 10- 5 (T r ..(i)- 1 R] / ?(tn) yr, on which 



a thin disk forms (jCannizzo. Lee fe Goodmanlll990l ; IMenou. Perna fe Hernquistll200ll ). Here 
T c6 is the typical temperature in the outermost disk annulus during this early phase, in 
units of 10 6 K, and Rd,s(to) is the initial disk radius in units of 10 8 cm. We set T c6 = 1 
throughout this paper. The subse quent evolution of the disk obeys the self-similar solution 
(jCannizzo. Lee fe Goodmanlll99Cll ). i.e., 



Md{t) 



M d (t ), 



< t < t Q , 
t>t , 



where the power index p = 1.25 ( IFrancischelli fe Wijersll2002l ). Md(t n) is a constant, which we 

normalize to the total mass of the disk Md = J °° M d dt, by M d / (5t ) (IChatterjee. Hernquist fe Narayan 



2000J). 



The subsequent evolution of the NS can be divided into three phases: 

1. The accretor phase - When the magnetosphere radius R m ~ 1.6 x l^fiB^ M l8 ^ cm 
is less than the corotation radius defined by R c = (GM/Q) 1 / 3 , the fallback matter is allowed 
to be accreted onto the surface of the NS. Here B 12 = B/10 12 G, M 18 = M/10 18 gs~\ and 
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f2 is the angular velocity of the NS. To evaluate the accretion torque, we divide this phase 
further into two sub-phases: a) If R m < the NS radius Rns, the torque exerted on the star 
is assumed to be J = MR]^ S Q K (R NS ), where Q%;(R) is the Keplerian velocity at R; b) if 
Rns < R m < Rci the accretion torque is given by J = 2MR^ n Q K (R m )[l — Q/Q K (R m )]. 



The accretion rate M of the NS is assumed to be limited by the Eddington accretion 
rate Me — 1 X 10 18 gs _1 for a 1.4 M & NS. However, during the early phase of accretion, the 
accretion rate may be so high that the radiation i s trapped in the flow and neutrino losses are 
important close to the NS surface (jColgatd Il97ll ; IChevalierlll989l ). The gravitational energy 
generated by infalling matter is mainly released by neutri n o loss es. Hence the Eddington 
limit does not work in this phase. According to IChevalierl (Il989l ). when the accretion rate 
drops to about 3 x 1O~ 4 M yr _1 , the reverse shock reaches the radiation trapping radius 
so that the photons can begin to diffuse out from the shocked envelope. A luminosity of 
Eddington limit is expected now, and this luminosity can reduce or even reverse the inflow 
outside the shocked envelope. Hence we adopt the Eddington limited accretion only when 
M < 3 x lO^Meyr- 1 . Otherwise we assume all the mass is accreted by the NS. 

Along with mass accretion, we assume that the magnetic field evolution follows the 
relation, B = Bq(1 + AM /Mo), where Bq is the initial magnetic field strength, and Mq is 
a parameter that determines the rate of decay with its typical valu e of ~ 10~ 5 — 10 -4 M(t 
( ITaam &: van den Heuvellll986l ; IShibazaki et al.lll989l ; lRomanilll990l ). 







2. The propeller phase - This phase begins when M decreases so that R c < R m . The 
infalling matter is assumed to be accelerat ed outward owing to the ce ntrifugal barrier, taking 
away the angular momentum of the NS (jlllarinov fc Sunyaevl 119751 ) . In this phase we use 
the same formula as in the accretor phase b) to estimate the propeller spin-down torque. 

3. The radio pulsarpha.se - As M decreases further, R m exceeds the light cylinder radius 
Ric = c/fl, or the disk becomes neutral, the accretion ceases and the NS becomes a radio 



pulsar. The neutral timescale of a fallback disc is t n 



[Rl 10 (t ) x 10 16 



L /M d (t )]~ 1/2 - 75 to 



( IMenou. Perna &: Hernquistll200ll ). where Rd,w is the disk radius in units of 10 10 cm. In this 
work, when either R m > Ri c or t > t n , we assume that the radio pulsar phase turns on and 
the subsequent NS spin evolution follows the magnetic dipole radiation prescription, until 
the NS crosses the "death line", i.e. B 12 /P 2 < 0.17. 

In Fig. 1 we show several examples of the calculated evolutions of the NS spin and 
magnetic field under typical conditions according the scheme described above. The model 
parameters are listed in Table 1. For the spin evolution, the thick solid, solid, and dotted 
lines represent the accretor, propeller, and radio pulsar phases, respectively. In case 1 we 
adopt the initial parameters as M d = 0.28 M , P = 14 ms, B = 3 x 10 12 G, M = 10" 5 M , 
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and Rd = 10 9 cm. The NS is found to be spun-up to a period of 2.88 ms in 72 yr, when 
the disk becomes neutra0. In the meanwhile, the magnetic field decreases dramatically to 
4.2 x 10 8 G. This is a strongly spun-up pulsar, demonstrating a possible way for the formation 
of PSR J1903+0327. In case 2 we lower the initial disk mass to O.O2M , while keeping all 
the other parameters same as in case 1. This leads to a lower Md(to), so that the spinning- 
up efficiency also becomes lower. The NS enters the radio pulsar phase with a relatively 
longer period of 11.26 ms and a higher magnetic field of 6.8 x 10 9 G. In case 3 we lower the 
initial disk radius Rd by a factor of 10 compared with in case 1. This gives a lower value 
of to and a higher value of Md(to), and the NS is spun-up to a period of 3.63 ms. In case 4 
we choose a slower decay of the magnetic field than in case 1 by setting a higher value of 
Mo = 10~ 4 M Q . It is shown that the period evolution remains almost unchanged, but the 
magnetic field decays to 4.2 x 10 9 G, about 10 times of that in case 1. A similar feature can 
be seen in case 5, if we increase the initial magnetic field to be B = 1 x 10 13 G. Finally in 
case 6 Md is further decreased to be 0.01 M Q , now the NS is spun-up to a period of 11.82ms 
during the first 39 yr and then spun-down during the propeller phase to a period of 11.92 ms 
during the subsequent 88 yr. After that time the disk becomes neutral and the radio pulsar 
phase begins. In the meanwhile the magnetic field decays to 6.6 x 10 10 G. This is a mildly 
spun-up pulsar. 



3. Population synthesis 



Based on the analysis in the above section, now we try to estimate the statistical proper- 
ties (e.g. the birth rate, total number, orbital period and eccentricity distributions) of RSR 
J1903+0327-like pulsars formed through the fallback disk accretion involved (FDAI) chan- 
nel in the Galaxy, by using an evolutionary population sy nthesis (EPS) method base d on 
the rapid binary stel l ar evo lution (BSE) code developed by iHurley. Pols &: Toutl (120001 ) and 
Hurley. Tout fc Pold (120021 ). This code incorporates the evolution of single stars with binary- 
star interactions, such as mass transfer, mass accretion, common-envelope (CE) evolution, 
SN kick, tidal friction, and angular momentum loss mechanisms (e.g. magnetic braking and 
gravitational radiation). 



We assume all stars are born in binary systems. The initial mass function (IMF) of 



x The neutral time is given by t n oc [R^ 3 M d (t )} 1/2 - 75 x R% 5 jMenou. Perna fc Hernquisfooil). T he value 
t n = 72 yr used here is much less than that (~ 10 3 ~ 10 4 yr) in iMenou. Perna fc Hcrnauist (2001 ). This is 
becau se (1) we calculated Md(to) by normalizing it to the total disk mass, while lMenou. Perna fc Hernquist 
(|200ll ) estimated it by assuming M d (to) = Md/t , whic h gives a larger value of M d ( to) th an ours; (2) the 
adopted value of Rd = 10 9 cm is larger than the value in lMenou. Perna fc Hernquistl ( 2001 ). 
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Kroupa. Tout fc Gilmord (119931 ) is adopted for the primary's mass (Mi) distribution. For 
the secondary stars (M 2 ) we assume a uniform distribution of the mass ratio M 2 /Mi between 
and 1. A uniform distribution of In a is adopted for the binary separation a. The star 
formation rate parameter is taken to be S = 7.6085 yr -1 , corresponding to a rate of ~ 
0.02 yr _1 for core-collapse SNe in our Galaxy, assume all the stars with masses > 8M die 
through SN explosions. 

It is known that a velocity kick is imparted to newborn NSs during core-collapse SNe, 
and a Maxwe llian SN kick distribut i on with dispersion V a ~ 200 — 400 kms -1 is usually 



adopted (e.g. Lyne &: Lorimer 1994 ; Hobbs et al.l I2OO5I ) . However, there is mounting 



evi- 



dence that some NSs may have received relatively small velocity kicks, of order less than 
50kms" 1 (jpfahl. Rappaport fc Podsiadlowskil I2OO2I ; [pfahl et alJliooa IPodsiadlowski et al. 



20041) . A probable mechanism th at causes small kicks is the electron capture (EC) SNe 
(IMiyaji et al.l Il98dl ; iNomotol Il984f). The explosion energies of such SNe are significantly 



lower than those of normal SNe dDessart et al.l 120061 ; iKitaura. Janka fc Hillebrandtl 12006 



Jiang. Chen fc Wangl 120071 ) . so that the nascent NSs receive a small kick velocity. It has 
been suggested that the ini tial masses of stars that explode via the EC SN mechanism 
may range in ~ 6 — 1 2 M G (IPodsiadlowski et al.l 120041 ; IKitaura. Janka fc Hillebrandtl 12006 



Poelarends et al.ll2008l ). Alternatively, core-collapse SN would be sub-energetic if the initial 
mass of its progenitor is greater than ~ 18 M . Simulations for a range of progenitor masses 
and different input physics suggest that the explosion energy decreases with increasing pro- 
genitor mass, whereas the binding energy of the star increases with increasing mass. This 
leads to a transition mass range of ~ 18 — 25 M abov e which the SN explosion is so energy- 
deficient that fallback accretion would be important (jFryerl 119991 ) . For the aforementioned 
reasons, we assume a Maxwellian kick distribution with V a = 20 kms -1 for the progenitor 
stars with initial masses of 8 < M < 12 M & and M > 18 M Q , otherwise we take V a = 265 
kms -1 . 

The fallback disk mass Md i s poorly known, but various arguments suggest tha t M<j < 
0.1 M may not be unreasonable (IChevalierlll989l ; lLin. Woosley &: Bodenheimerlll99ll ). There 
should exist a relation between the disk mass and the properties of the progenitors. Un- 
fortunately, there are big uncertainties in some key factors in the SN models, including the 
neutrino energy and transport, equation of state of gas at nuc l ear densities, convecti o n, ro- 
tatio n, neutrino emission asymmetry, etc (IMonchmeyerl Il99ll ; lYamada fc Satol Il995l ; iFryer 
1999|), each of which could strongly affect the SN energy and hence the mass of fallback. So 
we do not attempt to find a empirical formula to depict this kind of relation. Instead, we 
adopt simplified assumptions that sub-energetic SNe generally produce more massive fall- 
back disks, i.e., the logarithm of the disk mass log(Md/M ) is distributed uniformly between 
—5 and —1 for stars of masses 8 < M < 12 M or M > 18 Mq, and between —6 to —2 in 
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other cases. 

The initial magnetic fields are chosen so that log B is distributed normally with a mean 
of 12.4 and a standard deviation of 0.4. The magnetic field decay parameter M is set to 
be 10~ 5 — 10~ 4 M Q , and the disk radius Rd ~ 10 7 — 10 8 cm. The initial Pq is chosen to 
be distributed uniformly between 10 and 100 ms. In this work an MPSR is defined as a 
radio pulsar with spin period less than 10 ms, and we only take the FDAI MPSRs (either 
in binaries or isolation) into account. In the case of BMPSRs, if the companion star fills 
its Roche lobe before the pulsar passes the death line, we take the elapsed time before this 
moment as the upper limit for the lifetime of the pulsar. 

We divide the binary FDAI MPSRs into several categories according to their formation 
channels and the properties of the companion stars described as follows. The first case 
is that the pulsar is formed by core collapse of the primary star, and has an MS binary 
companion. Since the pulsar evolution can be strongly influenced by the stellar winds from 
the companion, which are intensive for massive stars, we only consider the MPSR-low-mass 
main sequence (MPSR-LMS) binaries, in which the companion star has a mass less than 
1.5 M Q . Alternatively the pulsar is formed by core collapse of the secondary, and already 
has an NS/BH binary companion. They are called NS/BH-MPSRs binaries. Finally, a single 
MPSR could form during either the first or the second SN that disrupts the binary. 

To investigate the influence of the model parameters on the final results, we design nine 
models by changing the values of the CE parameter a^J^I, the field decay parameter M , the 
lower and upper limit of the fallback disk mass Mdj ow , M^ up p in case of sub-energetic SNe, 
and the disk radius Rd (listed in Table 2). The calculated numbers and formation rates of 
various types of binary and single MPSRs in our Galaxy are listed in Table 3, and briefly 
described as follows. 

In model A (regarded as our control model) we adopt o>ce = 1-0, M = 10~ 5 M Q , 
M d = 10~ 5 — 0.1 M Q , and Rd = 10 8 cm. It is predicted that in the Galaxy there are ~ 200 
MPSR-LMS (PSR J1903+0327 may be one of them) and ten times more NS/BH-MPSR 
systems. Moreover, ~ 10 5 single MPSRs are p roduced. This number is comparable with 



estimated from observations (e.g. lLorimerlll995l ). suggesting that some fraction of the single 



MPSRs might form from the fallback disk accretion rather the traditional recycling channel 
in binary systems. Since their formation rate is 100 times lower than normal pulsars, it is 



2 The CE parameter acE is defined as the ratio of total binding energy of the envelope Ebi n d,i, and the 
difference between the initial and final orbital energies of the cores, i.e., acE = Ebind,i/ {E or bj — E or ^i), 
where E or b,f and E or b,i are the final and initial binding energy of the cores, respectively (see Hurley, Tout 
& Pols 2002 for details). 
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not surprising that no MPSRs have been detected in young SN remnants. 

In model B we increase the o>ce to 3.0, and find that higher o>ce produces more MPSR- 
LMS binaries. This is because binaries with a higher oice tend to survive during the CE 
evolution, especially for low-mass, close binaries which are more likely to be subject to 
coalescence during CE evolu tion. For the same reason, lower olce (model C) leads to fewer 



MPSRs (see iLiu fc Lill2006l . for more detailed discussions). In model D we adopt a slower 
magnetic field decay, resulting in a shorter duration of the accretor phase, hence significantly 
reducing the numbers of all three types of MPSRs. In model E we choose a lower value of 
Rd than in model A. As mentioned in §2, this will cause longer final periods of the pulsars 
(compare the curves of PI and P3 in Fig. 1), and produce fewer MPSRs. In model F the 
lower limit of the fallback disk masses is increased by a factor of 3, but the numbers of 
MPSRs change little, since MPSRs are produced only in the case of heavy disk accretion. 
This effect can also be illustrated in model G: when the upper limit of disk masses become 
half in model A, the numbers of MPSRs decrease by a factor of 10. To constrain the possible 
ranges of MPSR numbers, models H and I represent two extremely favorable and unfavorable 
cases for the production of MPSRs, respectively. It is found that the number of MPSR-LMS 
binaries may range from less than one to nearly one thousand in the Galaxy. 

We plot the eccentricity vs. orbital period distributions of the MPSR-LMS binaries for 
models A, B, and C in Figure 2. The distributions of the other three models are similar as 
that of model A. It is seen that in both models A and C most MPSR-LMS binaries are in 
wide orbits (P or b > 10 3 days) with eccentricities e > 0.3. For model B, the binary population 
seems to be distributed into two groups. One is similar as that in models A and C, the other 
stretches across large ranges of e (~ — 1) and P or b (~ 1 — 10 3 day) with a tendency of larger 
e accompanied with larger P orh . 



4. Discussion and Conclusions 

In this paper we suggest a born-fast formation channel for MPSRs in the Galaxy. In 
this scenario, a newborn NS may experience spinning-up and field decay phase via accretion 
from a fossil disk, established as a result of fallback following a SN explosion. For appro- 
priate choices of initial parameters, the NS could become a low field MPSR within < 10 3 
yr. In particular, this scenario can naturally explain the properties (millisecond spin and ec- 
centric orbit) of the recently discovered BMPSR PSR J1903+0327, without invoking stellar 
interactions within a GC or peculiar triple system. 

Our population synthesis calculations also suggest a population of single FDAI MPSRs 
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in the Galaxy. The formation of single MPSRs has not been well understood, especially how 
the pulsar has lost its binary companion after mass transfer. The conventional explanation is 
that the pulsar was recycled in a binary located in a GC, and then ejected out due to stellar 
encounters, or the high-energy radia tion from the pulsar has evaporated its companion (e.g. 
Ruderman. Shaham. fc Tayanilll989l) . The latter explanation may also require a GC origin 
of the binary (king et al.l l2005h . Although the fraction of BMPSRs and LMXBs in GCs 
is much larger with respect to the fraction in the Galactic field, it seems unlikely that all 
single MPSRs have originated from GCs. The fallback dis k accretion scenario migh t provide 
a possible way for the formation of not only single MSPs (IMiller &: Hamilton! 1200 ll ) but also 
planetary systems around them (ILin. Woosley & Bodenheimerlll99ll ). 



Obviously the formation scenario proposed in this work is subject to many uncertainties, 
so the results in Table 3 should be regarded as the optimistic cases. One of the biggest issues is 
to determine how much fallback material would have enough angular momentum at the time 
of collapse to allow the formation of a disk. Our population synthesis calculations suggest 
that there may be tens to hundreds of MPSR-LMS binaries lurking in our Galaxy. For PSR 
J1903+0327, a relatively large amount of fallback disk mass (> 0.1 M ) is required to produce 
the very short period (< 3 ms), implying that the progenitor star may have experienced an 
unusual SN-fallback history. It is interesting to note that t he large mass (~ 1 .74 M Q ) of the 
pulsar may suggest a massive (> 18 M Q ) progenitor star (jZhang et al.l 120081 ) and probably 
intensive fallback during the SN explosion. To form a centrifugally supported accretion disk 
around the NS, the specific angular momentum of the fallback matter should be larger than 
that needed for a circular orbit at the NS radius [GM^sRns) 1 ^ 2 ~ 1.4 x 10 16 cm 2 s _1 . Modern 
stellar evolution calculations suggest, however, that t he majority of massive star s lose a ngular 
momentum via magnetic processes very efficiently iHeger. Woosley fc Spruitl (120051 ) found 
that magnetic torques decrease the final rotation rate of the collapsing iron core by about 
a factor of 30 — 50 compared with the nonmagnetic counterparts. In their magnetic models 
there is a small amount of matter (~ 0.5 M for a 15 M Q star and ~ 1.0 M & for a 25 M & 
star) that has specific angular momentum greater than 10 16 cm 2 s _1 in the envelope of the 
pre-SN progenitor star. The realistic amount of fallback matter would be significantly less 
than these values because material with high specific angular momentum materials lies in 
the envelope of the progenitor star rather than in its inner core. Solutions of the problem 
of having enough rotation may involve tidal interactions in close binaries. If either star in 
the binary fills its Roche lobe, the spin-orbit coupling tidal torques would corotate the spins 
of two components with the orbital rotation, so that the progenitor star of the NS could 
be spun-up enough to allow the formation of the fallback disk. If this is true the realistic 
number of FDAI MPSRs should be less than our calculation, because not all the progenitors 
of the FDAI MPSRs have experienced the close binary evolution phases. 
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Other uncertainties include accretion from the SN fallback disk and magnetic field decay 
induced by accretion. The rate at which material falls back on the accretion disk depends on 
both the density profile of the star and on its angular velocity profile. Numerical calculations 
hav e been carried for a few mass and angular momentum configurations of massive stars 
(cf. IWoosley fc Blooml 120061 . and references therein), which generally suggest highly super- 
Eddington accretion with a short time. Still lack is a comprehensive picture of the conditions 
for fallback disk formation after SN explosions of massive stars. Highly super-Eddington 
accretion in a short period of ~ 100 yr in a radiation-trapped regime plays a key role in 
the formation of FDAI MPSRs. This distinguishes our model from the standard type of 
accretion expected in the recycling scenario, which lasts for ~ 10 s — 10 10 yr. In the standard 
pulsar recycling theory, the so called "spin- up line" in the B — P diagram shows the minimum 
period to which the NS can be spun-up in Eddington-limited accretion. This line is defined 
by the equilibrium period assuming the spin-up proceeding at the Eddington accretion rate. 
Although the pulsars in our scenario have experienced the highly super-Eddington accretion, 
since the mass accretion rate declines rapidly, the equilibrium period is never attained, unlike 
in the case of accreting NS in LMXBs. Thus our model does not necessarily predict any 
MPSRs sitting above the spin-up line. 

Nor is known how accretion-induced field decay occurs in NSs. One possible way is 
via ohmic dissipation of the accreting NSs crustal c urrents, due to the heating of the crust 



which in turn increases the re sistance in the crust (lRomanilll990l : iGeppert fc Urpinlll994 



Konar fc Bhaattacharyal 119971 ) . Alterna tive scenarios consist of screening or burying the 

magnetic field with the accr eted material ( Bisnovatyi-Kogan fc Komber gill 974 ; Taam Sz van den Heuvel 
19861 ; ICumming et al.ll200ll ). or outward moving vortices in the superfluid a nd super conduct 



ing c o re pushing magnetic fluxoids in to the crust during pulsar's spin-down (ISrinivasan et al. 



19901 ; iKonar fc Bhaattacharyal Il999l ) . In this work we require a fast field decay during the 
hyper-accretion phase (within < 10 3 — 10 4 yr). Screening or burying the magnetic field with 
the accreted material may be more appropriate for this scenario. 

Finally we move to the uncertainties is the CE evolution. We have used a constant CE 
parameter a CE to compute the orbital evolution during the spiral-in process, which is, how- 
ever, yery likely to change with the properties and evolutionary state of stars (llben &; Livio 
19931 ; iTaam fc Sandquistll2000l ). From Fig. 2 we find that model B with a C£ = 3 seems to be 
preferre d for the formati o n of PSR J1903. This value is be compatib le wit h other es t imate s 
by, e.g., 



van 



den Heuvell Jl994h iPortegies Zwart fc Yung' 
but in contradiction with Taam ( 19961 ). Sandquist et al. 



son ( 



1998 



998) and iKaloeeral (119 99). 



2000|), and lO'Shaughnessv 



(120081 ). Here the controversial issue is that there are no strict criteria for defining binding en- 
ergy of stellar envelopes and there is no clear understanding whe ther sources other than gravi- 
tational energy may contribute to unbinding common envelopes ( IHan. Podsiadlowski fc Eggleton 
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1995|). 

A distinct feature of the born-fast scenario is that the companion star of PSR J1903+0327 
is predicted to be a "young" (< 1 — 2 x 10 9 yr, i.e. within the characteristic age of the pulsar) 
MS star. Future optical, IR, and radio observations could present strong constraints on the 
nature of the optical counterpart, and thus verify or falsify the born-fast scenario. 

We thank X. -P. Zheng, H. -L. Dai and Z. -Y. Zuo for helpful discussions, and J. R. 
Hurley for providing the BSE code. It is grateful to the anonymous referee for his/her in- 
sightful comments that help improving the original manuscript. This work was supported by 
the Natural Science Foundation of China under grant numbers 10573010 and and 10221001. 
XWL thanks Huazhong Normal University for a postdoctoral fellowship. 
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Fig. I. — Period (bottom-left axes) and magnetic field (top-right axes) evolution of an NS 
with initial parameters list in Table 1. The curves labeled PI - P6 and Bl - B6 represent the 
evolutions of period and magnetic field, respectively. Numbers are suffixed to indicate the 
corresponding models. In the P — t plot, the thick solid lines, solid lines, and dotted lines 
represent the accretor, the propeller, and the radio pulsar phases, respectively. The magnetic 
field evolution is plotted until the radio pulsar phase begins. 
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Fig. 2. — Orbital period and eccentricity probability distribution of the MPSR-LMS binaries 
for model A, B, and C, from left to right, respectively. The distributions for the other six 
models (D-G, H, and I) are very similar to those for models A, B, and C. 
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Table 1: Initial model parameters for the NS-disk systems. 



Model No. M d (M ) R d:10 P (ms) B (10 12 G) M (1O" 5 M ) 



1 


0.28 


0.1 


14 


3 


1 


2 


0.02 


0.1 


14 


3 


1 


3 


0.28 


0.01 


14 


3 


1 


4 


0.28 


0.1 


14 


3 


10 


5 


0.28 


0.1 


14 


10 


1 


6 


0.01 


0.01 


14 


10 


1 



Model parameters for binary population synthesis. 


Model 


OtCE 


M (1O- 5 M ) 


Rd,8 


M dtlow (1O- 5 M ) 


M d)Upp (0.1 M ) 


A 


1 


1 


1 


1 


1 


B 


3 


1 


1 


1 


1 


C 


0.5 


1 


1 


1 


1 


D 


1 


10 


1 


1 


1 


E 


1 


1 


0.1 


1 


1 


F 


1 


1 


l 


3 


1 


G 


1 


1 


l 


1 


0.5 


H 


3 


1 


l 


3 


1 


I 


0.5 


10 


0.1 


1 


0.5 



-18- 



Table 3: Predicted numbers and formation rates in our Galaxy of various types of FDAI 
MPSR system^ 



"ft T 11 

Model 


item 


MPSR-LMS 


i\Tfi /-nTT i\ innn 

NS/BH-MPSR 


Single MPSR 


A 


number 


203 


3.3 x 10 3 


1.5 x 10 5 




rate 


6.0 x 10 


3.1 x 10 


3.1 x 10 


B 


number 


751 


3.7 x 10 3 


1.6 x 10 5 




rate 


1.3 x 10~ 6 


4.0 x 10~ 6 


3.2 x 10" 4 


C 


number 


192 


3.1 x 10 3 


1.5 x 10 5 




rate 


6.5 x 10~ 7 


3.0 x 10~ 6 


3.1 x 10^ 4 


D 


number 


8 


1.4 x 10 2 


1.2 x 10 3 




rate 


6.2 x 10~ 7 


3.3 x 10~ 6 


2.7 x 10~ 4 


E 


number 


83 


1.3 x 10 3 


4.9 x 10 4 




rate 


3.9 x 10- 7 


2.0 x 10~ 6 


1.9 x 10~ 4 


F 


number 


232 


3.7 x 10 3 


1.7 x 10 5 




rate 


6.9 x 10~ 7 


3.5 x 10~ 6 


3.4 x 10" 4 


G 


number 


28 


4.5 x 10 2 


1.2 x 10 4 




rate 


3.5 x KT 7 


1.7 x 10~ 6 


2.0 x 10" 4 


H 


number 


856 


4.2 x 10 3 


1.8 x 10 5 




rate 


1.5 x 10~ 6 


4.5 x 10~ 6 


3.6 x 10" 4 


I 


number 


0.3 


7 


38 




rate 


9.9 x KT 8 


4.6 x 10~ 7 


3.9 x 10~ 5 



